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I.   Introduction 

Theoretically,  the  oxidation  of  carbon  has  been  one  of  the  most 
intensely  considered  topics  in  the  field  of  combustion.   Experimentally, 
the  data  on  carbon  combustion,  especially  for  small  particles,  is 
lacking.  The  major  studies  and  experiments  with  carbon  sphere 
combustion  are  as  follows. 

In  1924,  Nusselt  proposed  his  "shrinking  drop"  model  for  the 
combustion  of  a  solid  sphere  such  as  wood  or  coal.  This  model  assumes 
that  oxidation  occurs  only  at  the  surface  producing  CO  and/or  CO.  and 
that  the  internal  temperature  history  is  governed  strictly  by  a 
1-dimensional  conduction  equation.  When  a  point  on  the  sphere  reaches  a 
characteristic  pyrolysis  temperature  it  is  assumed  to  have  burned. 
While  this  model  was  a  start,  it  failed  for  several  reasons.   It  ignored 
temperature  dependence  of  reaction  rates,  stoichiometry,  diffusion 
limits,  and  ash  effects. 

In  1931,  Burke  and  Schumann  assumed  that  pure  carbon  is  consumed 
due  to  the  reaction  with  C0„  at  the  surface  to  form  CO.  The  CO  thus 
formed  reacts  with  the  surrounding  oxygen  at  a  flame  sheet  to  form  C09. 
The  CO.  thus  formed  at  the  flame  sheet  then  diffused  back  to  the  surface 

to  sustain  the  process.   Burke  and  Schumann  showed  the  importance  of 
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stoichiometric  considerations  in  the  problem. 

In  1951,  Spalding  developed  the  principles  of  liquid  drop 

combustion  which  used  an  evaporative  vapor-phase  diffusion  flame  to 

describe  the  combustion  of  droplets.  While  this  model  adequately 

described  several  solid  and  liquid  combustion  systems,  it  failed  in  the 

case  of  carbon.  Due  to  the  high  sublimation  temperature  of  carbon  which 


is  3800°K,  it  really  does  not  have  a  significant  vapor  phase 
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concentration  at  realistic  combustion  temperatures. 

In  1958,  Coffin  developed  a  system  for  determining  the  steady  state 

burning  rate  of  solids  utilizing  mass  transfer  coefficients  rather  than 

vapor  phase  diffusion.   Using  this  method,  he  was  able  to  effectively 

model  for  the  cases  in  which  the  mass  transfer  processes  are 

4 
rate-determining . 

The  late  1970' s  brought  expensive  oil  and  a  greatly  renewed 
interest  in  coal  combustion.  This  factor,  coupled  with  an  explosive 
increase  in  the  availability  of  computational  facilities,  spawned 
several  computer  model  studies.  Principle  among  the  computer  model 
studies  was  the  one  by  Libby  and  Blake. 

Laurendeau,  in  his  tome  Heterogeneous  Kinetics  of  Coal  Char 
Gasification  and  Combustion,  covered  all  aspects  of  the  topic;  coal  and 
char  characteristics,  surface  mechanisms,  particle  reaction  models,  and 
kinetics  of  char  gasification  and  combustion.   Of  principle  interest  in 
this  study  is  the  comparison  of  the  Shrinking  Drop  Model  which  assumes  a 
strictly  surface  reaction  to  the  Progressive  Conversion  Model  which 
considers  pore  effects.   For  small  particles  on  the  order  of  tens  of 
microns,  there  is  considerable  question  as  to  which  is  the  proper 
mechanism. 

Libby  and  Blake,  using  a  computer  model  of  a  single  coal  particle 
that  considered  ash  effects  and  reaction  rates  as  well  as  diffusion, 
established  the  concept  of  an  extinction  diameter.  The  extinction 
diameter  is  the  particle  size  at  which  a  burning  particle's  heat  loss 
exceeds  its  heat  production.  At  the  extinction  diameter  the  particle 
extinguishes  itself.  ' 


Experimentally,  the  literature  is  sparse.   In  the  1930' s  studies 
were  done  by  Tu,  Hottel  and  Davis  to  determine  the  combustion  rate  of 
carbon.  Large  spheres  (2.5  cm  in  diameter)  and  carbon  samples  in  cups 
were  heated  in  furnaces,  the  mass  loss  rates  were  measured  and  the  gas 
profiles  around  the  carbon  bodies  were  sampled.  These  experiments 
validated  Burke  and  Schumann's  idea  of  the  flame  front  for  large 
particles.  However,  it  must  be  remembered  that  this  work  was  conducted 
with  large  particles.  '  ' 

Experimental  work  in  the  area  of  individual  coal  particle 
combustion  slowed  in  the  1940's  as  pulverized  coal  began  being  used  as  a 
fuel  in  boilers.   Interest  shifted  to  furnace  design  with  large  groups 
of  interacting  particles.  Sherman  in  1940  analyzed  the  space 
requirement  for  pulverized  fuel  combustion  using  an  entrained  flow  of 
coal  particles  through  a  furnace.    Even  though  data  was  needed  on 
individual  small  particles,  it  was  not  until  1976  that  work  in  this  area 
was  conducted. 

In  1976,  Ubhayakar  and  Williams  studied  the  burning  rate  of  small 
carbon  particles  (50  um  -  200  um) .  The  particles  were  suspended  on  a 
quartz  fiber  with  grease  and  ignited  with  a  pulse  from  a  ruby  laser.  As 
the  particle  ignited,  it  fell  from  the  quartz  fiber  and  was  followed  by 

a  high  speed  camera.  Additionally,  two  color  pyrometry  gave  a 

12 

temperature  history.  '   Results  showed  that  for  50-200  um  particles  the 

flame  sheet  idea  is  not  valid  and  the  flame  is  confined  to  the  particle 
surface.  Additionally,  particle  extinction  was  experimentally  verified. 

Since  the  study  by  Ubhayaker,  high  speed  cinematography  and 
two-color  pyrometry  have  been  used  in  several  studies  that  follow 
particles  through  a  furnace.   This  system  gives  burnout  times  and 
temperature  histories  for  50-100  um  particles. 
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Real  boilers  use  fuel  sieved  through  200  (74  Um)  or  400  (37  Um) 
meshes.  This  gives  a  size  distribution  weighted  in  the  20-50  um  range. 
Studies  have  not  been  conducted  on  this  size  carbon  particle  because  up 
to  now,  there  has  not  been  an  experimental  method  by  which  they  could  be 
measured.   In  the  past  15  years  a  new  version  of  an  old  device,  the 
modified  Milliken  Cell,  has  been  developed  that  can  isolate,  hold,  and 
weigh  these  sizes  of  particles. 

The  bihyperboidal  Milliken  cell  was  first  developed  in  1970  at 
Aberdeen  Proving  Ground  to  provide  a  method  for  isolating  and  observing 
20  to  100  um  aerosol  particles.    In  1981,  Arnold  developed  the  concept 
of  electron-stepping  as  a  method  of  mass  and  charge  measurement.  This 
concept  is  based  upon  the  fact  that  if  a  charged  particle  lost  an 
electron,  the  particle  mass  and  charge  can  be  determined  from  two 
voltage  measurements.  Arnold  applied  this  concept  to  a  flat  plate 
Milliken  cell  and  was  able  to  weigh  picogram  size  particles.    In  1982, 
Philip,  et  al.,  showed  that  the  technique  of  electron  stepping  could 
also  be  used  in  the  bi-hyperboidal  Milliken  cell  and  an  extremely 
sensitive  mass  isolation  and  measuring  device  was  available.    It  has 

the  capability  of  capturing  and  holding  a  particle  at  one  position  in 

—8       —1 2 
space  and  it  can  weigh  particles  in  a  mass  range  from  10  g  to  10   g. 

Since  1981,  the  bihyperboidal  Millikan  cell  or  quadrupole  has  been  a 

piece  of  equipment  in  search  of  a  use. 

The  objectives  of  this  study  were  to  extend  the  data  on  carbon 

particles  into  the  range  of  10  to  50  Um  and  to  study  the  absolute 

changes  in  these  particles  caused  by  the  combustion  process.  The 

Milliken  Cell  allows  a  direct  measurement  of  mass  and  size  changes  for 

these  particles.   Because  absolute  size  and  mass  change  data  are 


available  through  this  method,  evaluation  of  the  two  principal  small 
particle  combustion  models  is  possible.  Data  collected  in  this 
experiment  shows  that  for  20-50  urn  carbon  graphite  particles  a 
progressive  conversion  model  that  includes  pore  diffusion  is  preferred 
over  a  shrinking  core  model  that  ignores  pore  diffusion. 


II.   Experimental 

The  experimental  equipment  and  technique  for  mass  determination 
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followed  as  closely  as  possible  that  used  by  Philip,  et  al. 

A.   Equipment 

1.   Quadrupole  (Modified  Milliken  Cell) 

The  quadrupole,  shown  in  Figs.  1  and  2,  that  suspends  the  particles 
has  surfaces  defined  by  the  following  equations.  The  inner  surfaces  of 
the  endcap  electrodes  are  given  by 
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z2=ih  +  V*  (i) 

where  Z  =  4  mm  and  is  the  vertical  distance  from  the  geometric  center 

of  the  chamber  and  r  is  the  radial  distance  from  the  axis  of  symmetry. 

The  characteristic  dimension  of  the  chamber  is  Z  ,  half  the  distance 

o 

between  the  two  end-cap  electrodes,  measured  along  the  axis  of  symmetry. 
The  inner  surface  of  the  ring  electrode  is  described  by 

?        t-2 

z  2-i--V-  (» 

The  quadrupole  with  insulating  spacers  used  in  this  study  was 
fabricated  by  Temco  Corporation,  Danvers,  Mass. 

The  chamber  consisted  of  a  central  ring  electrode  of  aluminum  and 
two  teflon  rings  on  either  side.  The  two  aluminum  endcap  electrodes 
were  mounted  in  the  teflon  rings.  Ten  (1/8")  access  ports  were  drilled 
into  the  teflon  rings  between  the  ring  and  end-cap  electrodes.  All 
access  ports  could  be  sealed  by  flat  glass  windows  when  not  in  use.  All 
outside  metal  surfaces  on  the  quadrupole  were  covered  with  teflon 
shields  to  prevent  electrical  shock. 
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Fig.  I  Exploded  view  of  Quadrupole 
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Fig.  2    Cross-section   of   Ouodrupole 


2.  Circuitry 

The  circuitry  to  produce  the  ac  voltage  (with  dc  bias)  on  the  ring 
electrode  and  the  dc  voltage  across  the  endcaps  is  given  in  Fig.  3.  The 
capacitor-resistor  network  suppresses  induced  ac  voltages  on  the  endcap 
electrode  resulting  from  the  ac  voltage  on  the  ring  electrode.  The 
regulated  dc  power  supply  to  produce  V  across  the  end-caps  consists  of 
a  Pacific  Instruments  7104  PC  card.  This  card  is  capable  of  producing 
0-2000  V  across  the  end-caps  with  a  resolution  of  0.0005  V.  The  power 
transformers  were  Essex  Stancor  P-8151  2400  VAC  CRT  transformers.  The 
power  transformers  were  controlled  by  a  Powerstat  T/16  non-grounded 
variable  transformer.  The  ring  electrode  system  provided  an  rms  ac 
voltage  of  0  to  4800  V  at  60  HZ  with  a  resolution  of  about  50  V. 
Measurements  of  V  were  made  with  a  Keithly  192  digital  multimeter 
capable  of  reading  up  to  6J5  digits.  The  meter  was  protected  by  a  30  A 
fuse  on  the  negative  lead. 

3.  Illumination,  Particle  Sizing  and  Electron  Stepping 

Illumination  was  provided  by  a  Spectra  Physics  Model  102-3 
helium-neon  laser.  The  laser  was  mounted  on  a  single  axis  translational 
stage  and  attached  to  a  lead  brick.   This  mounting  coupled  with  a  beam 
steerer  allowed  the  laser  beam  to  be  directed  into  the  chamber  along  the 
asymptote  between  the  electrodes.  The  particles  were  observed  and  sized 
through  a  flat  glass  window  at  23°  to  the  illuminating  beam  using  a  Beck 
telescope  with  an  image  distance  of  200  mm,  a  50  mm  objective,  and  a  10X 
Beck  micrometer  eyepiece.   Sizing  was  accomplished  by  sweeping  the 
movable  reticle  across  the  particle  image.   Adjusting  the  ac  field  to 
align  one  edge  of  the  particle  with  the  stationary  stadia  line  gave  best 
results.  This  optical  system  allowed  a  resolution  of  2.5  ym  in  sizing 
particles. 
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The  UV  light  for  electron  stepping  was  provided  by  an  Oriel  6035 
low  pressure  mercury  calibration  lamp.  The  lamp  was  mounted  in  one  of 
the  ten  chamber  access  ports. 
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Fig.  3     Circuitry    for    Ouadrupole 
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4.  Combustion  Optics  (Fig.  4) 

Energy  is  delivered  to  the  particle  in  the  following  manner.  A  5  W 
Spectra-Physics  Argon  Laser  provided  a  beam  at  wavelength  of  514.5  nm. 
This  beam  was  split  and  focused  upon  the  particle  as  two  beams  180° 
apart  from  each  other  through  1/8"  holes  drilled  in  the  ring  electrode. 

The  beams  were  focused  by  100  mm  lenses  to  achieve,  for  a  0.5  watt  beam, 

8  6 

a  power  density  on  the  order  of  10  W/m2  and  a  heating  rate  of  10 

°K/sec.   See  Appendix  A.   Exposure  time  was  controlled  by  a  camera 

shutter.  The  shutter  was  calibrated  using  a  cadmium  photoresistor  and  a 

memory  oscilloscope.  Actual  times  are  shown  in  Table  1. 

Table  1.   Shutter  Calibration 


Shutter  Setting  Time  (actual) 

10  145.4  ms 

25  61.7  ms 

50  28.1  ms 

100  13.4  ms 

300  5.76  ms 


The  overall  experimental  schematic  is  shown  in  Fig.  5, 
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Fig.  4    Combustion    Optics 
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B.   Experimental  Procedures 

1,  Mass  Measurement 

a.)  Data  Collection 

The  ring  voltage  was  initially  adjusted  to  about  500  V.  The 
voltage  across  the  end-caps  was  set  to  its  minimum  value,  typically 
about  .025  V.  Particles  were  injected  into  the  chamber  with  a 
hypodermic  needle.  Due  to  the  applied  A.C.  field  in  the  ring  electrode, 
this  often  left  several  charged  particles  stably  oscillating  in  the 
quadrupole.  Then  a  particle  that  responded  suitably  to  the  balancing 
D.C.  voltage  was  balanced  at  the  center  of  the  chamber.  A  stationary 
particle  was  assumed  to  be  balanced  since  any  departure  from  the  center 
of  the  chamber  resulted  in  oscillations  of  the  particle.  The  ring-  ac 
voltage  was  then  lowered  for  a  few  seconds.  This  usually  allowed  the 
unbalanced  particles  to  fall  out  of  the  chamber  while  the  balanced 
particle  remained  near  the  chamber  center.  The  ac  voltage  was  then 
increased  again.  The  remaining  particle  was  more  carefully  balanced  by 
adjusting  V  and  then  lowering  the  ac  voltage  for  a  couple  of  seconds. 
If  the  particle  drifted  from  the  chamber  center,  it  was  clearly 
unbalanced.  The  ring  voltage  was  then  increased  again  and  V  further 
adjusted,  after  which  the  procedure  was  repeated,  until,  on  lowering  and 
raising  the  ac  voltage  the  particle  remained  stable.   The  UV  light  was 
then  turned  on  for  an  instant.   Electron  loss  was  verified  by  lowering 
the  ac  voltage  and  observing  particle  drift.  The  particle  was  then 
re-balanced  at  the  chamber  center  by  the  technique  already  described  and 
V  recorded  again.   This  procedure  is  repeated  a  number  of  times  and 
from  several  measurements  of  V  the  mass  and  charge  of  the  particle  can 
be  calculated. 
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b)  Mass  Calculation 

For  a  stationary  particle  at  the  center  of  the  quadrupole,  a  simple 

force  balance  describes  the  mass  of  the  particle  as 

C  q  V 
m  =  -j-g-i  (3) 

o 

where     q  is  particle  charge, 

C  is  the  geometric  constant  for  the  chamber  equal  to  0.4, 

g  is  the  gravitational  constant, 

Z  is  the  characteristic  dimension  of  the  quadrupole  Z  =  4  mm, 
and      V  is  the  potential  difference  across  the  end-caps. 

The  charge  on  the  particle  equals  the  number  of  excess  electrons 
times  the  charge  on  an  electron  so 

C  ne~V 

■--727*  •  (4) 

Collecting  constants 

m  =  K  n  V,  (5) 

where     K  =  %-  =  1.6356  x  10"15  g/V. 
gZo 

If  Vn  and  Vn+1  represent  the  balancing  voltages  for  a  particle 

containing  n  and  n+1  excess  electrons,  respectively,  then  by  solving  the 

previous  equation  at  each  balancing  voltage  the  particle  mass  is  given 

by 

V  V   , 
„   n  n+1  ... 

m  ■  K  — ~ —  •  (6) 

n    n+1 

-8 
For  relatively  large  particles  (10   g)  with  large  numbers  of 

excess  electrons  (10  )  over  a  small  range  of  voltage  (IV)  this  can  be 

approximated  as 
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y  2 

AV  for  a  single  electron 


The  problem  is  now  reduced  to  finding  the  voltage  change  for  a 
single  electron. 

Given  a  sample  series  of  voltage  measurements  a  first  order 
difference  table  is  set  up. 

312.824 

313.188  *  •;;: 

313.346   ™ 

313.463  *   -Li/  (8) 

313.559  >  -°J5 

313.601  >  ™ 

313.830  >  '  J 

20  21 
This  was  the  standard  method  used  by  Milliken  in  his  work.   ' 

From  here  Milliken  extracted  the  voltages  that  intuitively  looked  right. 

This  method  is  unsatisfactory.   It  is  too  subjective.  An  objective 

quantified  method  was  developed  to  deal  with  the  problem.  Possible 

values  of  AV  were  tested  for  fit. 

For  example,  assuming  various  values  of  AV  for  a  single  electron 

loss  for  the  given  data  and  dividing  the  differences  from  Eq.  (8) 

AV  .003  .004  .005  .006  .009 

98.00  73.50  58.80  49.00  32.67 

76.00  57.00  45.60  38.00  25.33 

39.00  28.25  23.40  19.50  13.00         (9) 

32.00  24.00  19.20  16.00  10.60 

14.00  10.50  8.40  7.00  4.67 

77.00  57.75  46.20  38.50  25.67 

Residual  .00  1.50  1.80  1.00  1.93 

From  these  tests  it  is  clear  that  the  voltage  change  for  a  loss  of 
a  single  electron  is  .003  volts,  since  that  value  gives  the  smallest 
residual. 
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In  this  case,  from  Eq.  (7) 

2,  ™„   ,  ,« 


m 


=  (313.391)  /. 003  x  1.635  x  10  (10) 


m  =  5.3546  x  10~8  g  .  (11) 


2.   Sample  Preparation 

The  carbon  particles  for  the  test  were  obtained  by  filing  a  piece 
of  CS-312  nuclear-grade  graphite  and  sieving  the  result  through  a  400 
mesh  (37  um) .  The  resulting  powder  ranged  from  10  to  50  um.  The 
density  of  the  particles  was  measured  as  follows.  A  101  mm  x  101  mm  x 
151  mm  block  of  graphite  was  weighted  at  2643.25  g.  The  resulting 
density  was  1.748  g/cm3.   It  was  assumed  that  the  ground  particles  have 
the  same  density  as  the  block.   Since  the  density  of  graphite  is  2.25 
g/cm3,  the  particles  were  about  22%  porous. 

Since  the  particles  were  porous,  an  attempt  was  made  to  backfill 
the  pores  with  oxygen.  Particles  were  placed  in  a  vacuum  for  two  days, 
then  under  a  pressure  of  two  atmospheres  in  oxygen  for  another  two  days. 
The  particles  were  placed  under  a  vacuum  or  "outgassed"  for  two  reasons. 

o 

First,  to  remove  all  gas  from  the  macro  (r  >  500  A)  meso  (20  <  r 

o  o 

<  500  A)  and  micro  (r  <  20  A)  pores;  and  second,  to  desorb  nitrogen  from 
the  active  sites  on  the  surface. 

After  out-gassing,  the  particles  were  soaked  in  oxygen  at 
2-atmospheres  to  refill  all  the  pores  with  oxygen  and  to  establish  a 
monolayer  of  oxygen  on  the  surface.  Once  treated,  the  particles 
retained  a  high  oxygen  concentration  for  several  hours  even  in  an  air 
ambient  for  the  following  reasons.   Oxygen  is  chemisorbed  and  tends  to 
stay  on  the  surface  more  readily  than  nitrogen,  which  is  only  physically 
adsorbed.  Additionally,  it  is  very  difficult  for  oxygen  to  escape  the 
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meso  and  micro  pores  in  the  absence  of  a  pressure  differential.   Since 
these  pores  account  for  approximately  50%  of  the  porosity,  oxygen  is 
slow  to  diffuse  from  the  particles. 

Ambient  conditions  in  the  quadrupole  were  changed  by  flowing  5 
cc/sec  of  gas  for  3  minutes.   This  accomplished  approximately  200 
changes  of  gas  inside  the  chamber. 

3.  Measurement  Sequence 

Measurements  were  made  in  the  following  sequence.  Particles  were 
injected  into  the  quadrupole.   An  ambient  atmosphere  was  established  in 
the  quadrupole  by  flowing  gas  through  it.   The  gas  was  turned  off.  A 
particle  was  brought  to  the  center  of  the  chamber  by  turning  off  the  AC 
voltage.  The  particle  was  stabilized  in  the  center  of  the  chamber  and 
it  was  sized  with  the  micrometer  eyepiece.  The  mass  of  the  particle  was 
determined  by  electron  stepping  as  described  in  Section  II-B-1.  A  laser 
pulse  was  impacted  upon  the  particle.  The  particle  was  sized.  The  new 
particle  mass  was  determined. 

C.  Uncertainty  Analysis 

Consider  the  sources  of  uncertainty  in  the  mass  calculation, 


m  =  K  V2/AV   , 

e 

«-  -  f  **■ +  ~  * +  m  m  •  (12) 


K  is  constant  for  all  mass  calculations  since  it  contains  only 
geometric  and  physical  constants.  Thus,  while  K  could  cause  a 
systematic  error  it  does  not  introduce  any  uncertainty  from  measurement 
to  measurement. 
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To  get  the  V  term,  V  is  measured  over  a  very  small  range  (1-2V)  at 
a  relatively  large  voltage  (300-500V) .  Therefore,  the  contribution  to 
the  uncertainty  of  measurement  is  at  most  around  1%.  For  example,  in 
the  worst  case  with  a  V  =  300  and  a  voltage  range  of  2V, 

—  5V  =  (K/AV  )  2V6V  .  (13) 

3V 

The  contribution  to  the  uncertainty  would  be  at  most  1.2%. 

The  major  source  of  uncertainty  is  the  determination  of  AV  . 
Typically,  values  of  V  range  from  .003V  to  ,030V.   Since  the  precision 
of  the  voltmeter  is  ±.001  V,  the  effects  of  a  1  mV  voltage  uncertainty 
can  be  large, 

|fv  5  AVe  "  (K^2)  W  6AVe  '  (14) 


If  AV  =  30  mV,  the  uncertainty  contribution  is  3.33%. 
If  AV  is  3  mV,  the  uncertainty  contribution  is  33.33%. 

The  degree  to  which  I  can  calculate  AV  is  the  key  to  the  accuracy 

of  the  measurement.   I  improve  my  certainty  of  AV  in  two  ways.  First,  I 

know  the  particle  size  ±  2.5  ym  and  I  know  the  density  of  the  material, 

this  gives  me  a  rough  estimate  of  the  mass.   Second,  I  rely  on  several 

measurements  of  V  in  the  determination  of  AV  .  The  rough  knowledge  of 

the  mass  prevents  gross  error  in  the  calculation.   The  use  of  several 

voltage  measurements  cuts  the  uncertainty  of  the  AV  to  .0005  V  and  the 

uncertainty  due  to  the  determination  of  AV  is  cut  in  half. 

e 

Depending  on  the  AV  magnitude  the  uncertainty  of  measurement  will 
range  from  ±  2%  to  18%  and  that  is  acceptable. 
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III.  Results 

A.  Power  Study 

The  first  series  of  runs  tested  the  effect  of  varying  the  energy 
delivered  to  the  particles.   Untreated  nuclear  grade  graphite  particles 
were  burned  in  air  with  decreasing  exposure  times  and  beam  powers, 
ranging  from  shutter  speed  015,  1.000  watts  to  shutter  speed  300,  0.007 
watts.  Three  significant  items  were  noted.  There  was  a  definite 
particle  swelling,  there  was  a  linear  relationship  between  initial  mass 
and  the  percent  mass  loss,  and  there  was  an  ignition  cut-off  at  a  low 
enough  laser  power.   Cut-off  power  was  shutter  speed  300,  0.015  watts. 

Particle  swelling  data  is  shown  in  Fig.  6.   Size  increase  is 
generally  about  20-30%  by  volume. 

The  linear  relationship  between  initial  mass  and  percent  mass  loss 
is  apparent  in  Fig.  7.   Additionally,  note  the  ignition  cut-off,  the 
particles  represented  by  E  and  F  neither  swelled  or  lost  mass.  Finally, 
the  particles  did  not  burn  to  completion,  they  reached  an  extinction 
point. 


Fig.  6  Power  Test  Swelling  Data 
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Fig.  7  Power  Test  Mass  Loss  Data 
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B.  Ambient  Effect 

After  the  power  test  was  completed,  the  effects  of  different 
ambient  atmospheres  were  explored.   Untreated  nuclear  grade  graphite 
particles  were  exposed  to  a  50  shutter  setting  laser  pulse  at  .50  watts 
in  ambients  of  pure  nitrogen,  50%  oxygen/50%  nitrogen,  and  pure  oxygen. 

The  swelling  data  is  presented  in  Fig.  8.   Trends  in  the  data  are 
inconclusive.  For  particles  that  ignited,  oxygen  concentration  seems  to 
have  little  effect  upon  swelling.   Oxygen  concentration  has  a  definite 
effect  upon  mass  loss  and  this  is  to  be  expected  from  purely  diffusional 
considerations.  Figure  9  shows  the  increase  in  mass  loss,  especially 
for  small  particles,  with  increasing  oxygen  concentration. 

Particles  in  pure  nitrogen  exhibited  neither  swelling  nor  mass 
loss,  one  could  say  they  did  not  ignite. 


Fig.  8  Ambient  Effect  Swelling  Data 
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Fig.  9  Ambient  Effect  Mass  Los  Data 
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C.  Soaking  Effect 

Swelling  of  pure  carbon  particles  was  unexpected.  While  swelling 

has  been  observed  in  bituminous  coal  particles  during  devolatilization, 

22 
it  has  not  been  noted  for  higher  grade  coals  or  graphite.    Swelling 

indicates  that  there  is  an  interstitial  or  pore  gas  effect.  This  effect 

was  tested  by  "soaking"  particles  in  oxygen.  Nuclear  grade  graphite 

particles  were  treated  by  exposing  them  to  a  vacuum  for  two  days  and 

then  back-filling  their  pores  with  oxygen  by  putting  the  particle  under 

2  atm  of  pure  oxygen  for  another  two  days.  Particles  were  then  exposed 

to  a  50  shutter  setting  laser  pulse  of  .50  watts  in  ambients  of  pure 

nitrogen,  air,  and  pure  oxygen. 

1.  Nitrogen  Ambient 

In  pure  nitrogen,  oxygen  soaking  had  no  effect.   Both  treated  and 
untreated  particles  failed  to  lose  mass  or  to  swell. 

2.  Air  Ambient 

In  air,  oxygen  soaking  had  a  significant  effect.  As  shown  in 

Figures  10  and  11,  swelling  increases  slightly  and  and  mass  loss 

significantly  increased  for  the  smaller  oxygen  treated  particles. 

3.  Oxygen  Ambient 

In  an  oxygen  ambient,  oxygen  soaking  has  little  effect  on  mass 
loss,  but  has  a  great  effect  upon  swelling.  Figures  12  and  13  show -the 
swelling  increase  and  the  mass  loss  comparison  for  the  treated  and 
untreated  particles. 


Fig.  10  Soaking  Effect  Swelling  Data  Air  Ambient 
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Fig.  11  Soaking  Effect  Mass  Loss  Data  Air  Ambient 


35 


E  2 

<  e 

Q  •-• 

U)  < 


o 
a 


a 


a 

03 


a 


a 

U3 


a 
in 


'SSOl   SSVH   lN33d3d 


a 


a 
m 


Fig.  12  Soaking  Effect  Swelling  Data  Oxygen  Ambient 


T 

,  — r 

r 

1 — 

■1 — 

— 1 — 

■  -1 1 

■ 

o 

» 

z 
o 
tn 

u 

z 

in 

• 

C\J 

0 

0 

0 

■ 

< 
a 

0 

• 

■ 

i- 
z 

0 

- 

■ 

UJ 

z 
a 

0 

- 

■ 

0 

0 

s 

-       0 

• 

0 

OXYCEN   TRE/ 
UNTREATED 

1 1— 

1 

i 

i 

1 

1 

1 

1                 1 

0) 


37 


.    GO 


.  in 


UJ 
N 

(/> 

-*  _J 

-•  < 


.    CO 


.    CM 


.  a 


V- 

t- 

s 

tn 

*t 

UJ 

§ 

< 

u 

z 

2 

3 

> 

O) 


3   g 


O) 


03 


to  in  ■«■  ro 

•3SV3a3NI    3ZIS 


(M 


Fig.  13  Soaking  Effect  Mass  Loss  Data  Oxygen  Ambient 
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IV.  Discussion 

The  Milliken  Cell  utilized  in  this  experiment  offered  a  unique  way 
to  examine  the  combustion  of  carbon.  Since  absolute  values  of  mass  and 
size  were  available,  direct  testing  of  models  is  possible.   In  the 
following  comparison  the  two  principle  models  for  small  particle  carbon 
combustion,  "the  shrinking  drop  model"  and  the  "progressive  conversion 
model"  are  compared  to  the  experimental  data. 

A.  Theory 

1.  Diffusional  Kinetics 

a)   Shrinking  Drop  Model 

Consider  the  laser  ignition  of  carbon  graphite.  The  sublimation 
temperature  for  graphite  is  3800  K  so  it  is  not  likely  that  a 
vapor-phase  diffusion  flame  can  exist  at  realistic  combustion 
temperatures  or  the  temperature  achieved  in  the  cell.   (App.  A).   Since 
there  is  not  a  vapor  phase  diffusion  flame  the  assumption  that  the 
chemical  kinetic  rates  are  much  faster  than  diffusional  rates  cannot  be 
made. 

Consider  a  carbon  surface  burning  in  a  concentration  of  oxygen  in 
the  free  stream  specified  as  C  .  The  burning  is  at  a  steady  mass  rate. 
If  the  surface  oxidation  rate  follows  first  order  kinetics  and  all 
reactions  occur  on  the  surface  then 

G   =  G./i  =  k  C   ,  (15) 

ox    f       s 

where  G   is  the  flux  of  oxygen  in  gm/sec  cm2,  G-  is  the  flux  of  fuel,  k 
ox  r 

is  the  heterogeneous  rate  constant  for  surface  oxidation  reflecting 
volume  to  surface  area  ratio;  and  i  is  the  mass  stoichiometric  index. 
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The  problem  is  that  C  ,  the  surface  oxygen  concentration,  is  unknown. 
But,  the  consumption  rate  of  oxygen  must  equal  the  diffusion  of  oxygen 
to  the  surface.  Thus,  if  h  is  designated  the  overall  convective  mass 
transfer  coefficient  one  can  write 


G   =  k  C  =  h   (C  -  C  )  .  (16) 

ox      s    o   °°    s 

Seeking  the  mass  burning  rate  in  terms  of  the  free  stream  oxygen 
concentration  C_.   It  follows  then  that 


kC     =h     C-h     C      ,  (17) 

S  O       °°  OS 


kC     +h     C     -  h     C      .  (18) 

S  OS  o       °° 


C     =   (h   /(k  +  h  )}  C^   ,  (19) 

so  o  °° 


or 


where 


and 


G       m   [tic  h  /(k  +  h  )}  CJ   -  K  Cm   ,  (20) 

OX  O  o  °°  °° 


K  =  k  h  /(k  +  h  )    ,  (21) 

O  o 


1/K  =   (k  +  h  )/k  h     =   1/h     +   1/k   .  (22) 

o  o  o 

When  the  kinetic  rates  are  large  compared  to  the  diffusion  rates 
K  =  h  ;  when  diffusion  rates  are  large  compared  to  the  kinetic  rate 

K  =  k.  This  describes  the  shrinking  drop  case  where  diffusion  into  the 

23 
pores  is  not  considered. 

Since  pore  diffusion  is  ignored  in  the  shrinking  drop  model  and  all 

reactions  are  assumed  to  occur  on  the  surface  the  following  results 

would  be  expected  in  terms  of  this  experiment.  For  graphite  with  a  low 

ash  content  the  particle  should  decrease  in  size  as  the  graphite  is 
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consumed  from  the  outside  in.  The  mass  loss  should  increase  with 
increasing  oxygen  concentration.  Treating  the  particle  by  back-filling 
its  pores  with  oxygen  should  have  no  significant  effect  upon  mass  loss 
or  swelling. 

b)   Progressive  Conversion 

Consider  the  situation  in  which  diffusion  to  the  particle  surface 
is  sufficiently  fast  that  it  is  not  the  controlling  rate.  For  the 
porous  medium,  carbon  is  being  consumed  within  the  pores  as  well  as  on 
the  surface.  The  surface  consumption  rates  are  therefore  controlled  by 
kinetic  rates;  however,  the  consumption  in  the  pores  are  controlled  by 
the  diffusion  of  oxygen  into  the  pore. 

Thus,  the  mass  consumption  rate  of  oxygen  in  terms  of  a  flux  of 
oxygen  must  be  that  which  is  consumed  at  the  surface  plus  that  which  is 
consumed  in  the  pores. 

G   =  k  C  +  D.  (3C./3n)   ,  (23) 

OX        S     1     1     s 

where  D.  is  the  internal  diffusion  coefficient,  C.  is  the  oxygen 
concentration  inside  the  particle,  and  n  is  the  normal  to  the  surface. 

Following  the  same  convention  used  in  the  shrinking  core 
development,  this  equation  can  be  written  in  the  form 

G   =  k  C  +  D.  (3C./3n)   =  k'  C   .  (24) 

OX       S     1     1     s        s 

Using  spherical  symmetry  the  internal  oxygen  diffusion  process  is 
described  by  the  equation 


D. 

i 


rd2C.   .  dC,. 
1  *  2   i 


dr     r  dr 


-  q  =  0  ,  (25) 
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where  q  is  the  oxygen  consumption  rate  per  unit  particle  volume.  With 
the  boundary  conditions, 

D.  4^-  =  h   (C  -  C  )  at  r  =  R  ,  (26) 

i  dr    o   °°    s 

fi£  =  o  at  r  =  0  .  (27) 

dr 

It  is  possible  to  express  the  quantity  q  as 

q  =  kS  C.  (28) 

where  S  is  the  internal  surface  area  in  a  particle  volume. 

The  solution  of  Eq.  (25)  in  terms  of  the  expression  as  given  in  Eq. 
(24)  results  in  the  following  expression  for  kf 

k'  =  k  +  X  D^^  [coth  (XR)  -  1/XR]  (29) 

where  R  is  the  particle  radius  and 

X  =  (S±  k/D^  (30) 

For  the  case  of  small  values  of  XR  (XR  <  0.55)  which  is  physically 

representative  of  a  small  particle  coth  (XR)  can  be  expanded  in  a  series 

24 

of  which  only  the  first  two  terms  are  considered  significant.    See 

App .  C . 

coth  (XR)  =  (1/XR  +  XR/3  -  -^1^-  +  ...)  (31) 

0  <  (XR)  <  ir 

Substitution  Eq.  (31)  into  Eq.  (29) 

kT  =  k  +  X2  D±  R/3  =  k[l  +  S±   R/3]  (32) 

The  expression  then  is  the  rate  constant  when  the  inner  pores 
participate.  Writing  this  in  terms  of  Eq.  (20). 
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G   =  {[k'h  /(k'  +  h  )}  C  =  K  C  (33) 

OX  O  o      »        oo 

K  =  k'h  /(k*  +  h  ).  (34) 

o       o 

In  terms  of  the  progressive  conversion  model  the  following 
experimental  results  are  expected.   Since  the  reactions  in  pores  are 
considered  and  the  particle  burns  throughout,  the  particle  should 
decrease  in  size  only  slightly  as  it  burns  from  the  inside  out.  The 
mass  loss  should  increase  with  increasing  oxygen  concentration. 
Treating  the  particle  by  back-filling  its  pores  with  oxygen  could  have 
an  effect  upon  mass  loss  or  swelling  by  initially  speeding  up  the 
reactions  in  the  pores. 

2.   Extinction 

Particle  extinction  in  the  sense  described  by  Ubhayakar  and 
Williams  and  Libby  and  Blake  describes  the  point  at  which  the  heat  loss 
rate  of  a  particle  exceeds  heat  production  in  the  particle  and  the 
particle  extinguishes  itself.   '    In  the  case  of  the  shrinking  drop 
model  which  is  the  model  used  by  the  above  investigators,  the  extinction 
condition  would  exhibit  itself  in  the  form  of  a  critical  diameter.  The 
critical  diameter  is  the  size  after  which  these  constant  density 
particle's  heat  balance  went  negative. 

The  progressive  conversion  model  would  differ  somewhat  in  its 
extinction  condition.   Since  the  particle  burns  from  the  inside  out  a 
critical  density  or  a  critical  mass  might  be  more  appropriate,  but  there 
would  still  be  a  point  at  which  the  heat  balance  went  negative  and  the 
particle  extinguished  itself. 
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B.   Comparison  with  Data 

1.  Power  Test 

The  power  test  mass-loss  data  is  consistent  with  either  of  these 
combustion  models.   Since  k  from  the  Arrhenius  Equation,  is  based  upon 
e-fold  temperature  dependence,  it  is  reasonable  to  see  an  ignition 
cut-off  in  either  model.  The  slope  of  the  data  also  seems  reasonable  as 
the  smaller  particles  are  closer  to  their  extinction  size  and  should 
therefore  have  less  mass  loss. 

The  power  test  swelling  data  is  not  consistent  with  shrinking  drop 
model.  According  to  the  Industrial  Engineering  Graphite  Handbook,  the 
CS  312  graphite  utilized  in  the  experiment  had  an  ash  content  of  0.09%, 

therefore  the  low  density  matrix  that  remained  after  the  laser  heating 

27 
was  essentially  pure  carbon,  not  ash.    With  a  surface  reaction  it  is 

difficult  to  imagine  any  mechanism  by  which  swelling  could  occur. 

The  power  test  swelling  data  can  be  explained  in  terms  of  the 

progressive  conversion  model.   Internal  pressure  developed  during 

combustion  in  the  pores  combined  with  the  decrease  in  structural 

integrity  as  the  particle  burns  could  result  in  swelling.   Additionally, 

local  extinction  conditions  could  be  reached  in  the  pores  as  the 

particle  density  decreased  giving  a  porous  mass  of  low  density  carbon. 

2.  Ambient  Test 

The  ambient  test  mass-loss  data  is  consistent  with  both  combustion 
models.  Mass  loss  increases  with  increasing  oxygen  concentration  and 
this  is  consistent  with  the  extinction  result  predicted  by  both  models. 

The  ambient  test  swelling  data  was  inconsistent  with  the  shrinking 
drop  model  and  could  be  explained  by  the  progressive  conversion  model 
for  reasons  previously  mentioned. 
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3.  Soaking  Test 

The  soaking  test  mass  loss  data  is  not  consistent  with  the 
shrinking  drop  model.   In  air,  the  mass  loss  increases  significantly  for 
treated  particles.   Since  the  shrinking  drop  model  ignores  the  pore 
effects  backfilling  those  pores  should  have  had  very  little  effect  upon 
mass  loss.   See  Appendix  D. 

The  soaking  test  results  can  be  explained  in  terms  of  the 
progressive  conversion  model.   Oxygen  reacting  in  the  pores  would 
further  open  the  pores  at  the  initiation  of  the  reaction.  The 
additional  internal  surface  area  would  then  increase  the  reaction  rate 
and  the  mass  loss. 

C.  Comparison  Result 

From  the  above  comparison,  the  progressive  conversion  model  is 
preferred.   It  can  be  used  to  explain  swelling,  and  it  can  explain  the 
increase  in  mass  loss  when  particles  are  treated  with  oxygen.  For 
carbon  graphite  particles  from  20  to  50  um  the  progressive  conversion 
model  that  considers  pore  diffusion  is  preferred. 
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V.  Conclusions 


The  bihyperboidal  Milliken  cell  has  been  established  as  a  viable 
apparatus  for  studying  heterogeneous  solid-gas  reactions  on  a  single 
particle  basis.   Using  this  device  data  on  carbon  particle  combustion 
and  extinction  was  extended  into  an  entirely  new  size  range. 

A  comparison  between  the  two  principle  models  for  small  carbon 
sphere  combustion  with  the  data  obtained  from  the  Milliken  cell  shows 
that  one  model  is  clearly  favored.  For  small  (20  ym  -  50  ym)  carbon 
graphite  particles  the  progressive  conversion  model  is  preferred. 

Finally,  a  standard  based  upon  a  known  substance,  CS312  Nuclear 
Grade  Graphite  with  known  chemical  composition  and  density  has  been 
established  for  future  char  combustion  work  with  other  carbon  bearing 
substances  such  as  coal,  wood,  etc. 

Two  further  studies  are  directly  suggested  by  this  experiment. 
First,  a  computer  model  study  using  a  constant  volume,  decreasing 
density,  progressive  conversion  approach  should  be  done.   Second, 
two-color  pyrometry  should  be  coupled  with  the  Milliken  Cell  to 
determine  particle  temperature  histories,  along  with  mass  and  size 
change  data.  This  additional  information  would  allow  a  more  complete 
elicitation  of  the  carbon  combustion  mechanism. 
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Appendix  A.  Particle  Heat-up  and  Steady-State  Temperature 

In  the  absence  of  a  chemical  reaction  the  steady  state  temperature 
can  be  approximated  using  a  simple  energy  balance. 

"Laser    conduction   ^radiation 

For  a  small  spherical  particle  in  a  large  enclosure  of  air,  the 
following  relationships  hold  at  the  particle  surface. 


4  tt  k  ,   (T  .     -  T  ,  .   J 
•       air   sphere ambient 

qcond  =        1/R  , 

sphere 

q   .  =  a4ir  R2  ,     e    .  .m      (T4  .     -  T   4) 
nrad        sphere  graphite    sphere    amb 

qT     =  it  R2  ,     (Beam  Power  Density)  (a    ..„  )  (Mirror  Loss) 
MLaser       sphere  J        graphite 

k    is  the  thermal  conductivity  of  the  air,  R  is  the  particle 

radius,  as  shown  in  Fig.  4,  the  focused  beam  waist  w  is  2.62  x  10   m. 

a 

Thus,  for  a  .5  watt  beam  the  power  density  is  2.32  x  10  w/m2.  The  beam 

reflects  from  4  mirrors  and  a  beam  splitter  on  the  way  to  the  particle. 

Assuming  a  10%  loss  at  each  surface,  the  mirror  loss  term  equals 

(.90)   =  .590.   Invoking  Kirchoff's  Law  and  assuming  absorptivity  equals 

emmissivity,  a  =  e  =  .92,  power  delivered  to  a  25  \im   particle  then 

becomes 

qT     =  tt(12.5  x  10"6)2(2.32  x  108)  (.92)  (.590) 
Laser 

_2 
=  6.1816  x  10   watts 


For  the  same  25  urn  particle 


q    ,  =  1.5708  x  10"5  k  .   (T  .     -  T  .) 
cond  air   sphere    amb 
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q   .  =  1.0242  x  10" 16  w/K4  (T  .    4  -  T   4) 
nrad  sphere     amb 


At  T  =  2000  K 


q    .  =  (1.5708  x  10"5)(137  x  10"3) (2000  -  300) 
cond 

-3 
=  3.66  x  10   watts 

qrad  =  1.0242  x  10~16  (20004  -  3004) 

-3 
=  1.6379  x  10   watts 


At  T  -  2500  K 


q     =  (1.5708  x  10"5)(222  x  10-3) (2500  -  300) 
cond 

-  7.6718  x  10"3  watts 

q   .  =  1.0242  x  10~16  (25004  -  3004) 
^rad 

=  4.0  x  10"  watts 


At  T  =  3000  K 


q     =  (1.5708  x  10"5)(466  x  10-3) (3000  -  300) 
cond 

=  1.9764  x  10   watts 

q   .  =  (1.0242  x  10"16)(30004  -  3004) 
rad 

_3 
=  8.29  x  10   watts 


At  T  =  3500  K 


q    ,  -  (1.5708  x  10"5)(1000  x  10~3)(3500  -  300) 
cond 

=  5.0266  x  10   watts 

qrad  =  (1'0242  x  10~16)(35004  -  3004) 
=  1.5369  x  10   watts 


From  this  very  rough  calculation  it  seems  that  the  temperature  will 
stabilize  at  a  maximum  temperature  of  3400°K.  However,  there  are  other 
loss  terms  that,  while  difficult  to  quantify,  are  definitely  present. 
The  surface  of  the  carbon  particles  are  not  smooth,  thus  both  the 
radiation  and  conduction  loss  terms  will  increase  due  to  increased 
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surface  area.   There  is  a  kinetic  energy  loss.  This  loss  manifests 
itself  in  the  oscillation  of  the  particle  about  the  field  null  point 
when  it  is  irradiated.   Oscillations  are  typically  500  Um  from  the  null 
point  at  60  HZ.  Finally,  there  is  the  question  of  aim,  these 
calculations  hold  only  for  the  center  of  a  Gaussian  Beam.  Considering 
the  additional  loss  terms  a  temperature  estimate  of  2000°K  to  2500°K  is 
probably  more  accurate. 

The  initial  heating  rate  can  be  calculated  from  the  energy  equation 

dT      qLaser 


dt   c„  m   *•  i 
P  partxcle 


6.1816  x  10~2 


(709)  1/6  ir(25.0  x  10_4)3(1.75  x  10-3  Kg) 


6.0897  x  106  °K/sec. 


At  this  rate  heat-up  would  be  very  fast  reaching  2000°K  in  0.33 


msec. 
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Appendix  B.   Ungraphed  Data 


Power  Test  Air  Ambient 


xlO"8 

xl08g 

xl0"8g 

Shutter 
Speed 

Power 

s 

0 

Sf 

M 

0 

Mf 

Am 

Am% 

15 

1.000 

18 

23 

3.37 

1.31 

2.06 

61.12 

15 

1.000 

12 

15 

2.42 

1.29 

1.13 

46.69 

15 

1.000 

16 

18 

4.60 

1.15 

3.45 

75.00 

15 

1.000 

10 

13 

2.26 

1.21 

1.05 

46.46 

15 

1.000 

14 

16 

3.85 

1.23 

2.62 

68.05 

15 

1.000 

12 

14 

3.48 

1.43 

2.05 

58.91 

25 

1.000 

13 

15 

1.47 

.87 

0.60 

40.82 

25 

1.000 

17 

21 

4.90 

1.12 

3.78 

77.14 

25 

1.000 

14 

16 

4.02 

1.00 

3.02 

75.12 

25 

1.000 

13 

17 

3.57 

1.14 

2.43 

68.06 

50 

1.000 

15 

18 

2.72 

1.37 

1.35 

49.63 

50 

1.000 

16 

19.5 

4.36 

1.23 

3.13 

71.78 

50 

1.000 

11 

14 

1.74 

1.00 

0.74 

42.52 

300 

1.000 

15 

17 

2.78 

1.30 

1.48 

53.23 

300 

1.000 

14 

16 

5.33 

.73 

4.60 

86.30 

300 

1.000 

13 

16 

2.83 

1.19 

1.64 

57.95 

300 

.500 

14 

16 

5.24 

.526 

4.72 

89.95 

300 

.250 

13 

15 

3.47 

1.03 

2.45 

70.41 

300 

.100 

15 

17.5 

3.01 

1.51 

1.50 

49.78 

300 

.027 

13 

15 

2.85 

1.54 

1.31 

46.06 

300 

.014 

15 

15 

5.80 

5.80 

0 

300 

.007 

13 

13 

2.82 

2.82 

0 

1  Size  unit  equals  2.5  Um 


Ambient  Test 
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Shutter 
Speed 


Power 


X108g 

S,      M 
f       o 


X108g 


X108g 
Am 


Am% 


PURE 

OXYGEN 

50 

.5 

12 

17 

1.14 

0.46 

0.68 

59.84 

50 

.5 

16 

18 

3.13 

0.85 

2.28 

72.95 

50 

.5 

12 

16 

1.14 

0.43 

0.71 

62.54 

50 

.5 

19 

23 

5.81 

0.73 

5.08 

87.42 

50 

.5 

16 

19 

1.87 

0.54 

1.13 

60.34 

50 

.5 

18 

21 

7.87 

1.01 

6.85 

87.10 

50%  OXYGEN /NITROGEN 


50 

11 

14 

2.84 

1.13 

1.71 

60.10 

50 

13 

16 

2.33 

1.03 

1.30 

55.82 

50 

13 

15 

3.16 

1.26 

1.90 

60.18 

50 

22 

25 

12.41 

1.38 

11.03 

88.87 

50 

14 

17 

4.06 

0.90 

3.16 

77.71 

50 

15 

17 
PURE 

2.35 

NITROGEN 

1.21 

1.14 

48.31 

50 

.5 

14 

14 

2.56 

2.56 

__ 

0 

50 

.5 

16 

16 

5.03 

5.03 

— 

0 

1  size  unit  equal  2.5  urn 
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Soaking  Test 


—8  —8        —8 

Shutter                        xlO  g  xlO  g  xlO  g 

Speed      Power     S      Sr      M  M.       Am       Am% 

r                  o      f       o  f 


PURE 

OXYGEN 

50 

.5 

11 

19 

1.10 

0.38 

.72 

65.4 

50 

.5 

16 

21 

3.62 

0.61 

3.09 

85.3 

50 

.5 

14 

20 

2.16 

0.72 

1.44 

66.7 

50 

.5 

10 

16 

1.09 

0.42 

0.67 

51.5 

50 

.5 

13 

18.5 

2.74 

0.86 

1.88 

68.6 

50 

.5 

15 

20 

4.70 

0.75 

3.95 

84.0 

AIR 


50 

.5 

11 

15 

1.13 

0.54 

0.59 

53.55 

50 

.5 

17 

21 

7.10 

0.56 

6.54 

92.17 

50 

.5 

16 

20 

7.06 

1.24 

5.82 

82.37 

50 

.5 

13 

16 

1.90 

0.69 

1.21 

63.91 

50 

.5 

14 

16 

3.00 

1.15 

1.85 

61.72 

50 

.5 

13 

16 
PURE 

3.57 

NITROGEN 

0.97 

2.60 

72.67 

50 

.5 

14 

14 

2.56 

2.56 

— 

0 

50 

.5 

13 

13 

1.62 

1.62 

— 

0 

50 

.5 

16 

16 

5.25 

5.25 

— 

0 

1  size  unit  equals  2.5  um 
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Appendix  C.  Approximation  of  XR 


X  =  [(Si  k)/D.]% 


S.  is  the  internal  surface  area  in  a  unit  particle  volume. 

i 

k  is  the  heterogeneous  reaction  rate  constant  for  surface  oxidation, 

D.  is  the  internal  diffusion  coefficient. 

1 

R  is  the  particle  radius. 


R 


Assume  a  25  urn  particle 

R  =  12.5  x  10~4  cm. 


si 


The  particle  has  .22  cm3  of  void  per  cm3  of  volume.  From  the 


Industrial  Graphite  Engineering  Handbook,  an  effective  pore  radius  of 

-4 
3.2  x  10   cm  can  be  assumed,  therefore  assuming  cylinders  winding 

28 
through  the  particle   where  L  is  the  total  length  of  the  cylinder, 


L  -  V     /i  ?, 
pores 


L  =  6.838  x  105  cm, 


and      S.  «  L  2*r  »  1.375  x  103  cm2/cm3. 

l 


Di 


From  the  Chemical  Engineering  Handbook,  pore  diffusion  for  a  gas 
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can  be  approximated  as 

-1 

/"  *  *  * 

D  -*    (1 


pore       t 

v4r 


f  ttM  "* 


2R6 
a 


% 


where 

X  is  the  internal  porosity  of  the  particles, 


55 


r  is  the  average  pore  radius, 

xis  the  tortuosity, 

6  is  the  absolute  temperature, 

R  is  the  gas  constant 

M  is  the  molecular  weight, 


and  Df  is  the  free  diffusion  constant. 


X-  -22 

r  -  3.2  x  10"4  cm  from  Table  5A.08.0^Qin  the  Industrial  Graphite 

Engineers  Handbook 

T-  3     Table  5A.08.02 

M  =  32  g/mole 

R  =  (8.314  x  107  erg/mol  K)(1.02  x  10"  g  £m2/erg) 

=  8.48  x  104  g  cm2/mol  k 


a  =  300  K  to  3000  K 
A 


Df  -  .14  cm2/S 


D 
pore    3 


22 


1  32       f  +  1  hi 

-4s  \TZ   ,a    lrt4N  ,onn.  f    .14  cmz  sec 


4(3.2x10  H)  12(8.48x10")  (300) 


D     =  D,  =  7.026  x  10~3  cm2/sec  at  300°K. 
pore    i 

Making  the  approximation  that  diffusivity  varies  with  the  square 
root  of  temperature, 


/T 

vv  -— w- 

/To 


D.(1500  K)  =  1.5711  x  10"  cm2/sec  . 
D.(3000  K)  =  2.2218  x  10-2  cm2/sec, 
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Using  the  kinetic  parameters  quoted  by  Libby,  the  intrinsic 
reaction  rate  can  be  calculated  from  the  following  expression 


R  =  (p  )  8710  g/cm2-sec-atm  exp 
e     o2 


-[18xl03/T] 


Assuming  1  atm  of  oxygen  for  simplicity, 
At  300 °K 


R  =  7.6269  x  10"23  g/cm2-sec, 


At  1500°K 


R  =  5.3516  x  10~2  g/cm2-sec, 


At  3000°K 


Now 


R  =  2.159  x  101  g/cm2-sec. 

R  =  k  C  and  k  =  R  /C 
e      s  e  s 

At  1  atm/300K  the  density  of  oxygen  is 

^-1 
'82.06  (300)' 


Pox  = 


MP 


RT 


Pox  =  1.299  x  10   g/cm3, 


(D(32)   J 


-3  _,__» 


-20 
and    k  =  5.8714  x  10    cm/sec 


At  1  atm/1500°K, 


Pox  =  2.699  x  10   g/cm3, 


and    k  =  1.9822  x  10  cm/sec. 


At  1  atm/3000°K, 


Pox  1.299  x  10   g/cm3, 


and    k  ■  1.662  x  10  cm/sec  , 
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Putting  terms  together 


300°K 


(1.375  x  103)(5.8714  x  10~2V 


(7.026  x  10'3) 


=  1.0637  x  10~7  cm-1 
XR  =  1.3296  x  10"10, 


1500°K 


(1.375  x  10"3) (2.699  x  10~4) 


1.5711  x  10 


-2 


=  4.8602  cm"1 
XR  =  6.075  x  10"3, 


3000°K 


(1.375  x  103)(1.662  x  105) 


h 


2.2218  x  10 


-2 


1.0142  x  105  cm-1 


XR  -  126. 


From  these  rough  calculations  pore  diffusion  should  be  significant 
through  at  least  2000°K  and  the  expansion  approximation  should  hold. 
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Appendix  D.  Mass  Loss  from  Pore  Oxygen 


-4 
Particle  diameter  =  25  x  10   cm. 

1        O  _Q 

Volume  -  4-  ir  D  -  8.1812  x  10   cm3, 
b 

Porosity  =  0.22. 
Volume  0-  =  Porosity  x  Volume 
=  1.7999  x  10-9  cm3. 


PV    (1  atm) (1.7999  x  10~9  cm3) 
n  = 


RT   (82.06  cnr'-atm/mole  K)  (300  K) 


n  =  7.3113  x  10"14  mole. 


2C  +  0  ■*■   2C0 


If  all  internal  oxygen  reacts 

-13 
1.4623  x  10    moles  of  carbon  will  react 


M   .     12.011  g/mole 
carbon 

Mass  of  carbon  reacting  equals 

-12 
1.7563  x  10    grams. 

—8 
Typical  total  mass  loss  equals  2.5  x  10   grams,  so  mass  loss  for 

pore  oxygen  would  be  undetectable  if  it  did  not  affect  the  combustion 

process  in  some  way  other  than  being  simply  consumed. 
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Appendix  E.  Arc  Carbon  Test 

Electrode  carbon  was  burned  in  air  to  provide  a  comparison  with 
Nuclear  Grade  Graphite. 

Density  Determination  of  Electrode  Carbon 
Using  a  Mettler  Balance 

m  =  0.2268  g 
Using  a  Micrometer 

Diameter  -  0.0984  in  =  .2499  cm 
Using  a  Meter  Stick 

Length  ■  3.50  cm 
Volume  =  0.1717  cm3 

p  =  Mass /Volume  =  1.33  g/cm3 

After  grinding  and  seiving  the  electrode  carbon  in  the  same  manner 
as  the  Nuclear  Grade  Graphite,  the  electrode  carbon  was  burned  in  the 
Milliken  Cell. 

Ranges  for  Nuclear  Grade  Graphite  particles  are  shown  by  lines. 
Swelling  data  for  the  electrode  carbon  showed  the  same  trends  as  the 
swelling  data  for  the  Nuclear  Grade  Graphite  with  an  increase  of  20-30% 
by  volume  (Fig.  14).  Mass  loss  for  the  electrode  carbon  tended  to  be 
less  than  the  mass  loss  for  the  Nuclear  Grade  Graphite  (Fig.  15). 
Electrode  carbon  mass  loss  ranged  from  25%  to  72%  versus  40%  to  90%  for 
the  nuclear  grade  graphite  particles  with  initial  masses  in  the  same 
range.  This  can  be  explained  by  the  progressive  conversion  model  in 
terms  of  the  difference  in  density  between  the  electrode  carbon  and  the 
nuclear  grade  graphite.   Because  the  electrode  carbon  is  less  dense, 
local  extinction  conditions  will  be  reached  sooner,  and  it  will  not  lose 
as  much  mass . 


Fig.  14  Electrode  Carbon  Swelling  Data 
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Fig.  15  Electrode  Carbon  Mass  Loss  Data 


63 


< 

U  K 

z 

Ul  Ul 

§  ffl 

a.  X 

t-  < 
u 

uj  < 


a 
o 


a 

0) 


a 

CD 


a 


a 

(D 


a 
in 


a 


a 
m 


SS01   SSVW   !N33U3d 


a 


64 


REFERENCES 

Husselt,  W.Z.,  Ver.  Deut.  Ing.,  68,  124  (1924). 

2 
Burke,  S.P.  and  Schumann,  T.E.W.,  Proceedings  3rd  International 

Conference  on  Bituminous  Coal,  2_,    485  (1931). 

3Spalding,  D.B.,  Fuel,  30,  121  (1951). 

4 
Coffin,  K.P.  and  Brokaw,  R.S.,  National  Advisory  Council  for 

Aeronautics,  Tech  Note  3929  (1957). 

Laurendeau,  Normand  M. ,  Progressive  Energy  Combustion  Science,  k_,    221 
(1978). 

Libby,  Paul  A.,  Combustion  and  Flame,  J38,  285  (1980). 

Libby,  Paul  A.  and  Blake,  Thomas,  R. ,  Combustion  and  Flame,  36,  139 
(1979). 

g 
Hottel,  H.C.  and  Davis,  H. ,  "Combustion  Rate  of  Carbon,  Ind.  Eng. 
Chem. ,  26,  889  (1934). 

9 
Hottel,  H.C.  and  Parker,  A.S.,  Ind.  Eng.  Chem.,  28.,  1334  (1936). 

Tu,  CM.,  Davis,  H.  and  Hottel,  H.C,  Ind.  Eng.  Chem.,  26,  749  (1934). 

Sherman,  R.A.,  Proc.  3rd  International  Conf.  on  Bituminous  Coal,  2_, 
370  (1932). 

12 

Ubhayakar,  Shivadev  R.  and  Williams,  Forman  A.,  Electrochem.  Soc, 

123,  No.  5,  747  (1976). 

13 

McLean,  W.J.,  Hardesty,  D.R.  and  Pohl,  J.H.,  Eighteenth  Symposium 

(International)  on  Combustion,  Pittsburgh:  The  Combustion  Institute, 

1239  (1981). 

14 

Seeker,  W.R.,  et  al.,  "The  Thermal  Decomposition  of  Pulverized  Coal 

Particles,"  Eighteenth  Symposium  (International)  on  Combustion, 

Pittsburgh:  The  Combustion  Institute,  1213  (1981). 

Frickel,  Robert  H.,  Shaffer,  Roy  E.  and  Stamatoff,  James  B.,  "Chambers 
for  the  Electrodynamic  Containment  of  Charged  Aerosol  Particles," 
Report  //ARCSI-TR-77041,  Chemical  Systems  Laboratory,  Aberdeen  Proving 
Ground,  Maryland  21010  (1970). 

Arnold,  Stephen,  Journal  of  Aerosol  Science,  _1°_,  49  (1979). 

Philip,  Mark  A.,  Gelbard,  Fred  and  Arnold,  Stephen,  Journal  of  Colloid 
and  Interface  Science,  41_,  No.  2,  507  (1983). 

1  Philip,  Mark  A.,  Ibid. 


65 
19, 


Philip,  Mark  A.,  ibid. 

Wlliken,  Robe: 
Chicago  (1918) 


20 

Milliken,  Robert  A.,  Electrons  (  +  and  -) ,  Chicago:  University  of 


21 

Milliken,  Robert  A.,  The  ELectron,  Chicago:  University  of  Chicago 


22M 


(1935). 
dean,  W.J.,  ibid, 


23 

Glassman,  Irvin,  Combustion,  New  York:  Academic  Press  (1977). 

Glassman,  Irvin,  ibid. 

25 

Ubhayaker,  Shivadev  R.,  ibid. 

26Libby,  ibid. 

27 

Industrial  Engineering  Graphite  Handbook,  National  Carbide  Division, 

Division  of  Union  Carbide  (1962. 

28 

Indus  t .  Eng .  Graph .  Hand . ,  ib  id . 

29 

Perry,  Robert  H.  and  Chilton,  Perry  C,  Chemical  Engineering  Handbook, 

5th  Edition,  New  York:  McGraw  Hill  Book  Company,  16-19  (1973). 

30 

Indust.  Eng.  Graph.  Hand.,  ibid. 


Single  Particle  Carbon  Combustion 


by 


Joseph  H.  Cerv 


B.S.,  United  States  Military  Academy,  1978 


AN  ABSTRACT  OF  A  MASTER'S  THESIS 

submitted  in  partial  fulfillment  of  the 
requirement  for  the  degree 

MASTER  OF  SCIENCE 


Department  of  Nuclear  Engineering 

KANSAS  STATE  UNIVERSITY 

Manhattan,  Kansas 

1985 


Abstract 

Experimental  data  on  the  mass  loss  and  size  change  of  a 
laser-ignited  particle  of  nuclear  grade  graphite  (20um  -  50pm)  in 
quiescent  mixtures  of  oxygen  and  nitrogen  at  room  temperature  are 
presented.  The  test  variable  was  the  oxygen  molar  fraction  of  the 
ambient  environment.  A  quadrupole  allowed  absolute  measurement  of  mass 
for  the  individual  particles,  and  a  linear  relationship  was  found 
relating  initial  mass  with  percent  mass  loss  during  the  reaction. 
Additionally,  size  measurement  showed  a  swelling  of  20-30%  during  the 
reaction.  Two  theoretical  models,  one  an  unreacted  shrinking  core  model 
assuming  a  receding  reaction  front  and  the  other  a  progressive 
conversion  model  that  accounts  for  particle  porosity  are  considered. 
Swelling  data  indicates  that  the  progressive  conversion  model  is  the 
preferred  one  for  (20  -  50  ym)  graphite  particles. 
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